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I.  INTRODUCTION 


An  analysis  of  the  performance  of  a  loeal/mobile  radio  communications 
system  utilizing  the  slotted  ALOHA  random  access  protocol  is  investigated  in  this 
thesis. 

A.  ALOHA  NETWORKS 

1.  Pure  Aloha 

In  Aloha  networks,  the  information  comes  in  packets.  In  a  conventional 
pure  ALOHA  network  the  users  transmit  their  packets  whenever  new  data  are  ready 
for  transmission.  When  two  or  more  users  transmit  their  packets  simultaneously,  a 
collision  will  occur  in  the  channel,  and  the  packets  may  overlap  partially  or 
completely  (Figure  1).  The  collision  destroys  the  information  in  the  packets,  and  the 
affected  users  have  to  retransmit  their  packets.  These  collisions  cause  a  major 
degradation  in  performance,  which  limits  the  maximum  channel  throughput  of  pure 
ALOHA  to  l/2e  or  0.184  [Ref.  1]. 

2.  Slotted  Aloha 

In  a  slotted  ALOHA  network,  time  is  segmented  into  slots.  Each  slot  has 
a  duration  equal  to  the  packet  transmission  time  or  length.  Users  are  required  to 
synchronize  the  transmission  of  their  packets  so  that  the  leading  edges  of  the 
packets  are  aligned  with  the  beginning  of  the  predetermined  time  slot  at  the  receiver 


Figure  1.  The  Collision  in  Pure  ALOHA  Network. 

Since  packets  arrival  at  the  receiver  are  synchronized,  when  more  than  one  packet 
is  transmitted  simultaneously,  a  collision  will  occur  in  the  channel,  and  the  packets 
will  overlap  completely  (Figure  2).  The  elimination  of  partial  packet  overlap 
substantially  improves  the  performance  of  the  network.  The  maximum  channel 
throughput  is  increased  to  1/e  or  0.368  (twice  the  maximum  channel  throughput  of 
pure  Aloha)  [Ref.  1].  This  maximum  channel  throughput  is  not  necessarily  the 


packet 


tagged  packet 


packet 


packet 


Figure  2.  The  Collision  in  Slotted  ALOHA  Network, 

throughput  attained  by  the  system.  In  heavy  traffic,  the  likelihood  of  collisions 
increases  significantly.  After  each  collision,  the  affected  users  have  to  retransmit  their 
collided  packets.  As  the  number  of  retransmitted  packets  becomes  large,  the  channel 
will  reach  a  point  where  a  collision  occurs  every  time  slot.  Hence,  no  packets  can 
be  transmitted  successfully  and  the  system  is  said  to  be  unstable. 
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B.  ALOHA  NETWORK  IN  LOGAL/MOBILE  COMMUNICATIONS 

The  Aloha  network  is  receiving  a  great  deal  of  interest  for  use;in  applications 
which  involve  large  physical  areas  (eg.  ships,  multi-story  buildings,  warehouses,  urban 
environments,  etc.)  and  require  random  multiple  access.  The  two  primary  limitations 
of  the  Aloha  network  are  its  poor  maximum  throughput  and  its  inherent  instability. 
However,  the  Aloha  protocol  is  still  desirable  because  it  is  easy  to  implement. 

C.  PACKET  CAPTURE  EFFECT 

1.  Advantages 

In  local/mobile  radio  communications,  each  transmitter  produces  a 
different  signal  power  level  at  the  receiver  due  to  the  differences  in  receiver- 
transmitter  range,  shadow  fading,  and  Rayleigh  fading.  The  arriving  packet  with  the 
highest  power  level  will  have  the  best  chance  to  capture  the  receiver.  As  a  result,  it 
is  possible  for  the  receiver  to  accurately  decode  one  of  the  packets  involved  in  a 
collision.  Since  the  collision  does  not  necessarily  destroy  all  involved  packets,  the 
channel  throughput  will  increase  considerably.  The  capture  effect  significantly 
improves  the  performance  of  the  Aloha  network  [Ref,  f  j. 

2.  Limitations 

The  advantage  of  improving  channel  throughput  is  achieved  in  favor  of 
users  near  the  base  station  and  at  the  expense  of  users  far  from  the  base  station.  The 
user  far  from  the  base  station  will  have  a  low  probability  of  capturing  the  base 
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station  receiver,  and  consequently  will  have  to  transmit  its  packet  more  often  than 
a  user  near  the  base  station  and  hence  suffer  a  long  packet  delay  [Ref.  2], 

D.  PERFORMANCE  MEASUREMENT 

Several  methods  have  been  considered  to  measure  the  performance  of  the 
Aloha  system.  They  are  packet  reception  probability,  the  correct  reception  of  the 
packet  address  header,  and  capture  measures  based  on  a  threshold.  Since  the  packet 
reception  probability  provides  the  best  measure,  this  method  is  selected  to  measure 
the  performance  of  the  system.  A  transmitted  packet  is  considered  to  have  captured 
the  base  station  receiver  only  when  the  entire  packet  is  received  correctly  in  the 
presence  of  one  or  more  interfering  packets  [Ref.  3]. 

E.  ORGANIZATION  OF  THE  THESIS 

The  analysis  is  divided  into  four  sections.  Section  II  investigates  the  capture 
effect  for  a  ring  model  in  local/mobile  communications.  Section  III  presents  the 
numerical  results  for  a  three-ring  network.  Section  IV  analyzes  the  system  throughput 
for  a  multiple  ring  Aloha  network  using  a  Markov  chain  model.  Finally,  discussions 
and  conclusions  are  presented  in  Section  V. 
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II.  PACKET  CAPTURE  PROBABILITY  OF  A  RING  MODEL 


The  analysis  of  the  packet  capture  effect  in  a  ring  ihodel  for  local  /mobile 
radio  systems  will  be  presented  in  this  section.  The  probability  of  correct  packet 
capture  will  be  evaluated  for  a  system  using  ideal  coherent  binary  phase-shift-keying 
(BPSK)  modulation. 

A.  SYSTEM  MODEL 

A  network  consisting  of  uniformly  distributed  users  in  ah  annular  ring  will 
be  studied.  A  uniformly  distributed  density  function  for  the  user  traffic  is 
assumed.  The  model  illustrated  in  Figure  3  is  the  system  model  for  this  study. 

The  traffic  generated  within  a  distance  r  of  the  base  station  is 

r 

Gir)  =  2n  jxg(x)  dx  (1) 

0 

where  r  is  the  distance  from  the  user  to  the  base  station.  The  user  traffic  density  is 

g(r)- 

The  total  offered  traffic  in  the  network  is  obtained  as  G  =  G(r  «>). 
Therefore 

•• 

G  =  2nj‘xg(x)dx  (2) 

0 
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Figure  3.  System  Mode) 

The  spatial  distribution  function  for  packet  generation  ( the  probability  of  a 
a  packet  generated  within  a  distance  R  <  r)  is  [Ref.  3] 


7 


r 

F^(r)  =  =  ^fxg(x)dx 


(3) 


Therefore  the  probability  density  function  for  packet  generation  is  found  to  be 


.  /V  2%  .. 

U  ir)-  -  -^rs(r) 


(4) 


The  probability  density  function  for  packet  generation  in  ring  R;,  defined  as 
ri<r<ri+i,  is 


Ur\R)  = 


Air) 


0  ,  otherwise. 


(5) 


Assuming  that  users  are  uniformly  distributed  in  an  annular  ring  of  inner  radius 
Rj  and  outer  radius  Rj,  we  have  the  user  traffic  den.uty 


8(r)  = 


(6) 


[  0,  otherwise 

Substituting  (6)  into  (4),  we  have  the  probability  density  function  for  packet 
generation  as  follows 


fR(r)  = 


2r 


0  ,  otherwise 


(7) 


The  distribution  of  packet  generation  of  a  ring  with  Rj  <  r  <  R2  is  shown  in 
Figure  4.  The  probability  of  a  packet  being  generated  within  a  ring  of  radius  R  <  rj 
is  given  by 


8 


Figure  4.  User  Traffic  Density. 


C-o 


which  can  be  evaluated  to  yield 


9 


w  = 


R^-Rl 


(9) 


Tlie  probability  of  a  packet  generated  within  a  ring  of  radius  R  <  is 


r^-rI 


(10) 


Substituting  (10)  into  (5)  and  (7),  we  get  the  probability  der^city  function  for 
packet  generation  in  ring  Rj 


(11) 


where  i  equals  0,  1,  2,...,n-l  and  Rj  =  {q  ^  r  ^  Ti+i} 


B.  BPSK  BIT  ERROR  PROBABILITY 

We  assume  that  users  communicate  with  the  base  station  using  coherent  binary 
phase  shift  keying  (BPSK).  For  ideal  coherent  binary  phase  shift  keying  (BPSK) 
modulation  in  a  Rayleigh  fading  channel,  the  received  signal  of  a  particular  user  0 
in  the  presence  of  n-1  interfering  users  is  given  bv  ^assuming  the  background 
Gaussian  noise  i-,  much  smaller  than  the  interference) 

foit)  =  ejsj^t)  +  Vc.  e~\(t)  (12) 
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where  Sj(t)  are  BPSK  signals,  the  Cj’s  are  Rayleigh  distributed  signahamplitudes,  and 
the  6i’s  are  the  uniformly  distributed  phases  of  the  interfering  signals  (we  set  ©q  = 
0  and  assume  that  6j,  i=l,  2,...,  n-l  are  relative  to  this  reference). 

At  the  sampling  time  the  decision  variable  of  BPSK  receiver  is 

ii-i 

''o  =  ^0^0  +  fl,.  cos6.  (13) 

where  a;  =  {±1}  and  are  independent  and  identical  distributed  random  variables 
(i.i.d.).  The  second  term  on  the  right  S'de  of  ♦he  above  expression  is  the  sum  of 
independent  Gaussian  random  variables. 

The  received  signal  power  is 

S  =  E  ((Co  a/}  =  E  (cq}  (14) 

The  received  noise  power  is 

N  =  E  {(v^-  Co  fl/}  (15) 

Substituting  equation  (13)  into  equation  (15),  we  have  the  received  noise  power 
as  follows 


H-l 


cos0,.)^) 

/-I 


(16) 


Noting  that  the  mean  of  cos^©;  is  1/2,  we  find 
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(17) 


)7=:lE£fcf> 

^  i-1 


Hence,  the  received  signal  to  noise  ratio  is 


SNR  =  -  = 


S  2E 


N 

Y:E{cf) 

i-1 


(18) 


The  bit  error  probability  of  user  0  in  the  presence  of  n-l  interfering  users  for 
BPSK  in  a  Rayleigh  fading  channel  is  [Ref.  4] 


Pu  =  - 


1  - 


1 


1+ 


\  SNR  J 


(19) 


The  near/far  effect  is  modeled  by  assuming  that  the  received  power  w  conforms 
to  an  inverse  a  power  propagation  law;  that  is 


w  =  cr-*  (20) 

where  r  is  the  distance  of  the  user  from  the  base  station  and  a  is  the  propagation 
constant.  Generally,  3<a<5  for  ultra  high  frequency  (UHF)  applications.  The 
constant  c  depends  on  the  gains,  heights,  and  power  of  the  transmitter  and  receiver 
antennas.  Without  loss  of  generality,  we  assume  that  c  is  normalized  to  unity  for  all 
users;  that  is,  users  are  taken  to  be  identical  and  use  omnidirectional  anteimas  with 
radiation  pattern  maxima  in  the  horizontal  plane. 

From  (20),  the  received  power  of  user  i  is 
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Wj  =  E  icf)  =  Tj"* 


Substituting  (21)  into  (18),  we  have 


SNR  = 


2w„  2rn 


n-l  »-l 

E*,  Eo* 


H-l  /  f\-lt 


E 

i-l 


V^oy 


and  the  bit  error  probability  is 


Pu  =  - 


f,-  >  1 

n-l 

i-E 

i-l 

/O; 

-a 

which  can  be  written  as 


Pu  = 


1  - 


«-l  /_  \a 


\  1.1 


For  numerical  analysis  purposes,  we  select  the  propagation  constant  a 
Hence, 


1 


Pb  =  - 
^  2 


1  - 


1 


\  /-l 


The  bit  error  probability  for  the  case  of  M-level  diversity  is  [Ref.  3] 


(26) 


k-O 


'M-UkYl+\i 

.  k  ll 


\i 

> 


where  the  value  of  /x  is  given  by 

(27) 

C.  CODING  FOR  PACKET  CAPTURE 

Since  all  packets  are  transmitted  with  equal  power  and  because  of  the  neaf/far 
effect  and  the  Rayleigh  fading  channel,  the  operating  signal-to*noise  ratio  can  be  very 
low.  To  overcome  this  problem,  we  propose  a  combination  of  diversity  and  forward 
error  correction.  The  ^’ite^bi  decoding  scheme  is  employed.  This  scheme  proposes 
maximum  likelihood  decoding,  which  is  relatively  easy  to  implement  for  codes  with 
small  memory  orders.  A  binary  convolutional  code  with  rate  R  =  1/2  and  constraint 
length  7,  and  Viterbi  hard  decision  decoding  is  considered.  Define  the  input 
modulation  as  M-ary  and  the  output  modulation  as  Q*ary.  When  binary  codiiig  is 
used,  the  modulator  has  only  binary  input  (M  =  2).  Similarly,  when  binary 
demodulator  output  quantization  is  used  (Q  =  2),  the  decoder  has  only  binary 
output.  In  this  case,  the  demodulator  is  said  to  make  hard  decisions  (M  =  Q).  To 
simplify  the  implementation,  most  coded  digital  communication  systems  use  binary 
coding  with  hard-decision  decoding.  Since  errors  at  the  output  of  the  Viterbi  decoder 
occur  in  burst  even  if  the  errors  at  the  input  of  the  decoder  are  independent,  the 


n-l 

A 

1-E 

_2 

j-l 

/ij 
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probability  of  correct  packet  reception  cannot  be  evaluated  by  assuming  independent 
errors.  It  has  been  shown  that  probability  of  correct  packet  reception  is  lower 
bounded  by  [Ref.  3] 

Pc  ^  (1-P/  (28) 

where  L  is  the  packet  length  and  Pu  is  the  union  bound  first  event  error  probability 
[Ref.  5].  For  hard  decision  decoding,  a  tight  upper  bound  on  Pg  is  given  by 


d=df 


(29) 


where  a^  is  the  number  of  code  words  of  Hamming  weight  d,  d^  is  the  free  distance 
of  the  convolutional  code,  and  P,j  is  the  first  event  error  probability  of  a  code  word 
of  weight  d.  The  first  event  error  is  made  with  probability  [Ref.  5] 


P4‘ 


«»(a+l)/2 


,  d  odd 


U  d\  jjdp. 

2U/2> 


e=UdlZ 


d\ 


(30) 


Pm  ^  E  \  Pm  H-Pm^'^  »  ^ 


The  bound  of  (29)  can  be  further  simplified  by  noting  that  for  d  odd 


e=(.d*m\^) 


(31) 


Expressions  involving  P^j^j  can  be  moved  out  of  the  summation  since  P^m  is 
independent  of  the  summation  index  e.  As  a  result,  the  upper  bound  for  P^  is 
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p,  <  a-Ptu)^  E  r  (32) 

e-(<f+lV2 

Taking  the  lower  limit  of  e  at  0  instead  of  (d+l)/2,  we  obtain 

p,  <  pSS  e(‘']  (33) 

«»o  \ 

The  sum  can  now  be  evaluated  to  obtain  an  upper  bound  on  of 

P,  <  2'*  pS  (34) 

Similarly,  for  d  even  an  upper  bound  on  is  found  to  be 

^  *  E  ff].n«a-)’»«)®  (33) 

This  expression  can  be  simplified  in  a  manner  similar  to  that  used  to  obtain 
(34)  to  get  an  upper  for  P^  when  d  is  even 

P,  <  pZ  (X-PiJP  (36) 

which  is  identical  to  (34). 

The  correct  packet  capture  probability  q(n)  is  obtain  by  averaging  the 

r 

probability  of  correct  packet  reception  over  all  possible  location  of  n  users.  Thus 

^(")  =  / /•••/ (37) 

where  fR(r)  is  given  by  (7).  Now  substituting  (30)  into  (42),  we  obtain  a  lower  bound 
for  the  correct  packet  capture  probability  as 


(38) 


in.  NUMERICAL  RESULTS 


A  numerical  analysis  has  been  performed.  Programs  using  Matlab  provide 
numerical  results  for  the  above  correct  packet  capture  probability  for  user  0,  Eight 
interfering  packets  have  been  used  in  the  analysis.  As  the  number  of  the  interfering 
packets  increases  the  computation  time  becomes  very  lengthy.  It  took  approximately 
24  hours  to  compute  q(l),  q(2)...q(6).  Therefore,  the  simplified  P<jin  (36)  is  used  to 
calculate  q(l),  q(2),...  q(7),  q(8).  The  result  for  q(n)  using  the  simplified  bound  of  P^j 
in  (36)  yields  results  that  are  veiy  close  to  those  obtained  with  (30).  For  each 
additional  user  after  eight  interferers,  the  results  do  not  change  significantly.  The 
source  codes  in  Appendix  A  are  executed  to  obtained  the  results  which  are  listed  in 
Table  II,  III,  IV  ,  and  IV.  The  results  and  source  codes  using  (30)  are  listed  in 
Appendbc  B  for  reference  purposes. 

Numerical  results  are  obtained  for  the  specific  case  of  L  =  400  bits  and  a  = 
4.  The  levels  of  diversity  considered  are  M  =  1,  2,  and  4.  A  convolutional  code  of 
rate  1/2,  constraint  length  7  is  used.  The  values  of  a^  and  d  are  listed  in  Table  I. 
[Ref.  3] 
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Table  I.  THE  COEFFICIENTS 

OFaa  ANDd 


d 

11 

10 

38 

12 

193 

14 

1,331 

16 

7,275 

18 

40,406 

20 

We  will  study  a  ring  model  for  the  following  cases. 

A.  A  ONE-RING  NETWORK 

A  one-ring  is  used  for  the  network.  The  user  Tq  locates  in  a  range  from  = 
0  km  to  R2  =  50  km,  any  other  users  will  vary  in  the  same  range  with  user  rQ.  The 
probability  of  correct  packet  capture  for  one-ring  is 

5050  50 

q(n)  =  //.../(I  -?//,('■(»'•, (39) 

0  0  0 

where 
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(40) 


50^ 

0 


,  ^50 

,  otherwise 


where  j  =  0,  1,  2,  3, n-1.  The  calculated  values  for  q(n)  are  listed  in  Table  II. 


Table  11.  SIMULATION  RESULTS  FOR  TOE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  q(n)  FOR  ONE-RING 
NETWORK 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

^9 

q(V); 

q(8) 

M  =  1 

1 

B!!M 

.1116 

.0587 

.0382 

.0334 

II 

1 

.3676 

.2142 

.1116 

.0741 

.0649 

M=4 

1 

.6318 

.3964 

.2801 

.2132 

09 

B.  A  THREE-RING  NETWORK 

The  following  is  an  analysis  of  a  three-ring  model.  Ring  one  is  {0  km  <  r  <  25 
km};  ring  two  is  {25  km  <  r  <  38.5  km}  and  ring  three  is  {38.5  km  <  r  <  50  km}. 

1.  Ring  one 

The  probability  of  correct  packet  capture  for  ring  one  is 

50  50  25 

0  0  0 

where 


20 


fMR,)  = 


(42) 


—  ,  0^  r.  i  25 
25^  ® 

[  0  ,  otherwise 


,  0  ^  r,  s  50 
502  ^ 

0  ,  otherwise 

V  ' 


(43) 


where  j  =  1,  2,  3  The  calculated  values  of  q(ri)  are  listed  in  Table  III. 

Table  IIL  SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  RING  ONE 


q(l) 

IM 

q(3) 

.q(4) 

q(5) 

iIBB 

in 

q(8) 

M=1 

1 

.6758 

.4412 

.3221 

.2532 

.2080 

.1770 

.1543 

M=2 

1 

.8579 

.6960 

.5722 

.4801 

KB8W 

.3396 

.3141 

M=4 

1 

.9042 

.8099 

.7197 

.6350 

.4901 

.4615 

2.  Ring  two 

The  probability  of  correct  packet  capture  for  ring  two  is 

50  50  38.5 

g(n)  =  f...f  /(l-P//;j(rjR2)4(ri)../^(r„.i)dro..Jr„.i  (44) 

0  0  25 


where 


fjiW  = 


2ro 

38.52-252 

0 


,  25  s  To  S38.5 
,  otherwise 


(45) 
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(46) 


n)  = 


50^ 

0 


0^  rj  i  50 
, otherwise 


where  j  =  1,  2,  3,...n-l.  The  calculated  values  of  q(n)  are  listed  in  Table  IV. 


Table  rv.  SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  RING  TWO 


q(l) 

q(2) 

q(3) 

q(4) 

m 

q(6) 

WM 

q(8) 

M=1 

1 

0 

0 

0 

DBI 

0^ 

M=2 

1 

1!B!W 

.0019 

0 

OB 

0 

1 

.4731 

.2581 

.1080 

.0446 

.0063 

.005 

3.  Ring  Three 

The  probability  of  correct  packet  capture  for  ring  three  is 

50  50  50 

«W  =  /■■•/  /  (47) 

0  03tJS 

r 

where 


/rW  = 


]  502-38.5^ 
0 


38.5  s  5  50 
,  otherwise 


(48) 
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j  — 4  *  0  ^  0  ^ 

fR(^P  =  1  50"  ^ 


[  0  ,  otherwise 

where  j  =  0,  1,  2,...n-l.  The  calculated  values  of  q(n)  are  listed  in  Table  V. 


Table  V. 


SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  RING  THREE 


q(l) 

wm 

q(3) 

q(4) 

q(5) 

ilMIM 

q(8) 

B81 

1 

0 

0 

QHH 

0 

0 

0 

M=2 

1 

.0133 

nm 

0 

0 

0 

0 

0 

fl 

1 

.2302 

.0210 

.0008 

0 

0 

0 

0 

Plots  of  q(n)  versus  n  for  case  1)  and  each  of  the  rings  for  case  2)  are  shown 
in  Figures  5,  6,  7,  and  8. 
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Figure  5.  Probability  Of  Correct  Packet  Capture  For  One-Ring  Network 
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Total  Number  of  Packet 
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uer  of  Packet 


Figure  7.  Probability  Of  Correct  Packet  Capture  For  Ring  Two 
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Figure  8.  Probability  Of  Correct  Packet  Capture  For  Ring  Three 
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IV.  SYSTEM  THROUGHPUT  FOR  A  THREE  RING  ALOHA  NETWORK 


The  analysis  of  a  ring  model  for  mobile/local  radio  communications  with 
correct  packet  capture  probability  has  been  performed  in  the  previous  section.  Tl  e 
numerical  results  show  that  users  who  are  located  close  to  the  base  station  are  more 
likely  to  capture  the  base  station.  The  throughput  for  each  ring  will  be  calculated  by 
using  a  Markov  chain  model.  The  Markov  model  provides  a  detailed  model  of  the 
network’s  dynamic  behavior;  therefore,  it  is  a  desirable  model  for  the  analysis. 

The  probability  q;  that  one  of  i  colliding  packets  will  survive  a  collision  and  that 
the  other  i-1  packets  will  need  retransmission  later  has  been  calculated  in  Section  III. 
The  previous  system  model  of  an  Aloha  network  with  N  terminals  (users)  sharing  one 
channel  for  transmission  to  one  receiver  (base  station)  is  used.  It  is  assumed  that 
positive  acknowledgements  from  the  receiver  are  never  lost,  so  the  user  will  not 
retransmit  unnecessarily.  This  assumption  has  been  shown  to  be  very  realistic  since 
modern  mobile  networks  usually  have  to  guarantee  a  ni  ety  five  percent  coverage  of 
the  serving  area.  The  values  of  qj  are  calculated  in  Section  III,  and  it  is  assumed  ♦i.at 
qj  does  not  equal  to  zero  for  i  greater  than  one.  This  assumption  also  depends  on  the 
specific  properties  of  the  network.  In  an  undisturbed  channel,  qj  is  equal  to  1  and  qo 
is  equal  to  0  [Ref.  8]. 
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During  any  time  slot  each  terminal  (user)  will  be  in  any  of  three  states:  the 
original  state  O,  the  transmission  state  T,  or  the  retransmission  state  R.  The  terminal 
will  be  in  state  O  when  there  is  no  packet  ready  to  transmit  or  if  a  new  packet  has 
just  arrived  while  it  is  waiting  for  the  next  time  slot.  The  terminal  will  be  in  state  T 
when  it  is  actually  transmitting  a  packet,  successfully  or  not.  The  terminal  will  be  in 
state  R  after  it  has  transmitted  a  packet  unsuccessfully  or  when  it  is  waiting  for  a 
future  time  slot  for  retransmission.  The  terminal  will  be  designated  as  an  0-,  T-,  or 
R-terminal  depending  on  which  state  it  is  in.  It  is  assumed  that  the  process  which 
generates  new  packets  will  be  halted  in  a  terminal  while  it  has  a  packet  to  transmit. 
The  possible  state  transitions  for  an  individual  terminal  are  shown  in  Figure  9. 

It  is  assumed  that  changes  of  state  can  occur  at  the  beginning  or  at  the  end  of 
the  time  slot  in  the  Aloha  network.  An  0-terminal  can  either  remain  in  state  O  (do 
nothing)  or  change  into  a  T-terminal  at  the  beginning  of  a  time  slot.  Upon  changing 
into  a  T-terminal,  it  starts  the  transmission  of  a  new  packet.  The  generation  rate 
which  a  terminal  changes  from  state  O  to  state  T  is  denoted  Po.  An  R-terminal  can 
either  stay  in  state  R  (do  nothing)  or  change  into  a  T-terminal  at  the  beginning  of 
the  time  slot.  In  the  latter  case,  it  starts  the  retransmission.  The  retransmission  rate 
when  a  terminal  changes  from  R  to  T  will  be  denoted  p^.  At  the  end  of  a  time  slot, 
a  T-terminal  can  either  go  to  O  after  a  successful  transmission  or  to  R  after  an 
unsuccessful  transmission. 
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A  slotted  -  Joha  network  can  be  described  as  a  homogenous  finite  Markov 
process  with  N+1  possible  states  (0,  1,  ...N)  corresponding  to  the  number  of  R- 
terminals  at  the  end  of  each  time  slot. 

Given  0  <  n  <  N  and  0  <  m  <  N,  let  denote  the  state  transition  probability 
that  the  network  will  move  from  state  n  to  m  in  one  step.  There  are  three  different 
types  of  transitions: 

Type  1,  in  which  the  state  decreases  more  than  one.  This  is  not  possible 
because  only  one  packet  will  be  transmitted  successfully  in  one  time  slot. 

Type  2,  in  which  the  state  decreases  by  one  with  probability  qj.  This  will  happen 
if  and  only  if  there  is  no  0-terminal  and  i  R-terminals  change  into  state  T  for  1  <  i 
<n  (one  packet  is  successful). 
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Type  3,  the  state  moves  to  m  >  n.  This  will  happen  with  probability  qi+n,.n+i  if 
and  only  if  m-n+ 1  0-terminals  and  i  R-terminals  change  into  state  T  for  0  <  i  <  n; 
or  with  probability  l-qj+m-n  if  and  only  if  m-n  0-terminals  and  i  R-terminals  change 
into  state  T  for  0  <  i  <  n  (no  packet  successful). 

It  is  assumed  that  the  terminals  are  operating  independently.  The  generation 
rate  and  the  retransmission  rate  are  denoted  po  and  p^  respectively.  The  probability 
of  correct  packet  capture,  qj,  was  obtained  in  Section  III.  The  value  of  is 
calculated  as  [Ref.  8] 

0  ,  if  m  <  n-1 

‘‘  Sf”]  •  W-  (1-P/'^-  Q  ,  if  m  =  n-1 

i*l  \  * 

a-pf--'.  p, 

<-o  \  O 

I  \  /I  \ 

To  calculate  the  equilibrium  state  occupation  probabilities  an  arbitrary 
positive  number  is  assigned  to  tt^*.  The  is  calculated  as  follows  [Ref.  6] 


(50) 

,  if  m  n 


-  /  n-1  ' 

«  1  *  * 

'^n  • 

^n.n-1  \  «=0 


Normalizing  the  tt^*,  one  gets 


/«o 


(52) 
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The  expected  throughput  of  the  network  in  state  n  (the  probability  to  transmit 
a  packet  successfully  in  a  time  slot)  is  found  to  be 

jb-0  \iy-n-K)  1-0  V  *  / 

The  throughput  for  a  ring  is 

(50 

n-0 

The  system  throughput  is  shown  in  Figure  10.  The  curves  SI,  S2,  and  S3  are  the 
throughputs  of  ring  one,  two,  and  three,  respectively.  The  curve  S1+S2+S3  is  the 
overall  throughput  of  the  network.  The  source  code  for  Figure  10  is  listed  in 
Appendix  C. 


V.  CONCLUSIONS 


At  first,  an  introduction  to  Aloha  networks  have  been  made.  The  slotted  Aloha 
network  is  selected  in  this  study,  since  its  system  throughput  is  double  the  pure  Aloha 
system  throughput.  Then,  the  probability  of  correct  packet  capture  of  slotted  Aloha 
is  studied  for  mobile/local  radio  communications.  An  application  of  slotted  Aloha 
in  a  ring  model  has  been  developed.  In  this  model,  the  collisions  do  not  necessarily 
destroy  all  the  involved  packets.  With  the  help  of  this  model,  the  probability  of 
correct  packet  capture  has  been  investigated.  It  has  been  shown  that  for  a  ring 
network  the  system  throughput  is  greater  than  1/e,  and  the  ring  network  is  more 
stable  under  overload.  The  numerical  results  in  Section  III  have  been  plotted  in 
Figures  5,  6,  7,  and  8.  The  probability  of  correct  packet  capture  of  one-ring  and 
three-ring  networks  have  been  studied.  The  numerical  results  and  its  plot  in  Figure 
5  show  that  as  the  number  of  users  grows  from  one  to  six  the  probability  of  correct 
packet  capture  drops  rapidly.  For  each  additional  user  c/ler  sbc  users,  the  probability 
of  correct  packet  capture  does  not  decrease  as  quickly  as  with  only  a  few  users. 
When  the  number  of  users  is  large,  the  system  reaches  a  steady  state.  The  probability 
of  correct  packet  capture  increases  as  M  (diversity  level)  increases. 

Numerical  results  for  a  three-ring  network  show  that  when  the  users  are  located 
in  the  inner  ring  (e.g.  ring  one)  they  will  have  a  better  chance  to  capture  the  receiver. 


34 


As  the  users  move  to  the  outer  ring  (eg.  ring  three),  the  probability  that  they  will 
capture  the  base  station  decreases  rapidly.  The  system  throughputs  for  each  ring  have 
been  calculated  in  Section  IV  based  on  the  probability  of  correct  packet  capture 
obtained  previously. 
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APPENDIX  A 


As  mentioned  in  Section  III,  equation  (30)  is  used  to  calculate  q(n)  in  Tables 
II,  III,  IV  ,  and  V.  The  source  codes  for  the  results  shown  in  Tables  H,  HI,  IV,  and 
V  are  listed  below.  Since  q(6),  q(7),  and  q(8)  require  a  long  computation  time,  the 
program  is  written  to  calculate  each  value  of  q(n)  separately  for  large  value  of  n  (n 
>  6). 

The  source  code  for  the  numerical  results  of  Table  II  is  listed  below.  The 

number  of  users  (nuser)  is  initialized  to  2, 3,...,  and  8  to  calculate  q(2),  q(3) . q(8), 

respectively.  In  the  listed  program  the  maximum  number  of  users  (nuser)  is  equal  to 
eight.  For  a  one  ring  network,  Rj  and  R2  are  varied  from  0  km  to  50  km. 
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%(’  this  program  calculates  q(8)  for  one  ring  nehyork  *) 

%(*  the  coefficient  values  for  d  and  a^  *) 
d  =  [Id  12  14  16  18  10]; 
ad=[ll  38  193  1331  7275  40406]; 

the  number  of  interval  for  integration  *) 
interv=5; 

■%(^  the  values  for  Rj  and  Rj  *) 

=rl=[0  0]; 
r2=[50  50]; 

%(*  set  value,  of  a,  L,  number  of  user,  arid  ml  for -M  *), 
alpha =4; 

L=400; 

nuser=8; 

ml=[12:4j; 

%(*  calculate  the  grid  points  of  integration  *) 

deit'a(i) = (r2(l)-rl(l))/interv; 

delta(2) = (ri(2)-rl(2))/ihterv; 

r(i,l)=rl(l)+.delta(l)/2; 

r(2,,l)=ri(2)+delta(2V2; 

for  i=.2:iiiteiv 

r(l,i) = r(l,i-l) + deUa(l); 
end 

for  i=2:intefv 

r(2^i)=r(2,i-l)+delta(2); 

end 

%(*  calculate  the  probability  density  function  f^  equation  (7)  *) 
for  i=l:interv 

fr(l,i)=2*r(l,i)/(r2(l)"2-rl(ir2); 

end 

for  i=l:inteiv 

'  fr(2,i) =2*r(2,i)/(r2(l)"2-rl(2r2); 
end 


p(*  this  loop  calculate  q(8)  for  all 
%value  of  M=l,,2,  and  ,4,  equation  (38)  *). 
for  12-1:3 
ih=inl(i2); 
sum=0; 
for  a=l;inte:!^ 
for  b=l:mtefv 
for  c=l:inten' 
for  f=l:interv 
for  g;=l;irite.rv 
for  fi=l:iriteirv 
for  k=l:interv 
for  o=l:ihterv 

index=ia  b  c  f  g  h  k  oj; 
muy=(); 

%(*  calculate  value  for  equation'  (27)  *) 
for  il=l:nuser-l 

muy=mu}  +  (r(l,index(l)/f(2,iridex(il:+ l)))'^alpha; 
end 

muy = i/sqrt(l + muy); 

%(•  calculate  value  for  Pb^  epation  (25)  *) 
pb=(i'-muy)/2 

%(*  calculate  the  value  for  Pt,M>  equation' (26)  *) 

pbni=0; 

for  cnt=0:m-l 

pbm = pbm + (comb((m-l + crit),cnt)  *  (((1 + inuy)/2)'^cnt)); 
end 

pbm = (pb^'m)  *pbiri; 

%(•  calculate  the  value  for  for  d  even,  equqatioh  (36)  *) 
pd=  (2*sqft(pbm*(l-pbm))).^d; 

%{*  calculate  P^.,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,pe=l;,end 
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%(*  calculate  the  'integration  for  ■q(h);  'equati6h'(38)r*): 
ihc=(l-pe)';L*fr(i,o)*ff(2,a)*frtiB)*fr(2,c)*fr(^^^ 
inc=inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum = Slim +inc; 
end 
end 
end 
end 
end 
end 
end 
end 

sum=sum*delta(2)"'(nuser-l)*delta(l); 

end 


The  source  code  for  the  numerical  results  of  Table  III  is  listed  below.  For  ring 
one,  Rj  varies  from  0  km  to  25  km,  and  R2  varies  from  0  km  to  50  km.  This  source 

code  is  used  to  calculate  q(8).  To  compute  q(2),  q(3) . q(n),  the  number  of  users 

(nuser)  is  equal  to  2,  3,...,  n,  respectively. 


) 


%{*  this  program  calculates  q(8)  for  ring  one  *) 

%(*  the  coefficient  values  for  d  and  a^  *) 
d  =[10  12  14  16  1840]; 
ad =[11  38  193  1331  7275  40406]; 

%{*  the  number  of  interval  for  integration  *) 
interv=5; 

%{^  the  values  for  Rj  and  R2  *) 
rl=[0  0]; 
i^=[25  50]; 

%{*  set  value  of  o:,  L,  number  of  user,  and  ml  for  M  *) 
alpha =4; 

L=400; 
nuser=8; 
ml=[12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l) = (r2(l)-f  1  ( l))/intery; 

delta(2) = (r2(2)-rl(2))/interv; 

r(l,l)=rl(l)+delta(l)/2; 

r(2,l)=rl(2)+delta(2)/2; 

for  i  =2:interv 

r(l,i) =r(l,i-l) + delta(l); 
end 

for  i=2:interv 

r(2,i)  =r(2,i-l) + delta(2); 
end 

%{*  calculate  the  probability  density  function  f^  eq.  '  .ion  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(l)-2); 

end 

for  i=l;interv 

fr(2.i) =2*r(2,i)/(r2(ir2-rl(2n); 
end 


40 


%(*  this  loop  calculates  q(8)  for  all 
%value  of  M=l,  2,  and  4,  equation  (38)  *) 
fori2=l;3 
m=ml(i2); 
sum=0; 
for  a=  l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:interv 
for  h  =  l;interv 
for  k=l:interv 
for  o=l:interv 

mdex=[a  b  c  f  g  h  k  o]; 
muy=0; 

%(*  calculate  value  for  /i,  equation  (27)  *) 
for  il=l:nuser-l 

muy = muy  +  (r(l,index(l)  /  r(2,index(il + l)))'^alpha; 
end 

muy  =  1  /sqrt(l + muy); 

%(*  calculate  value  for  Pi,,  equation  (25)  *) 
pb=(l-muy)/2 

%(*  calculate  the  value  for  equation  (26)  *) 

pbm=0; 

for  cnt=0:m-l 

pbm = pbm + (comb((m-l + cnt),cnt)  *  ((( 1 + muy)/2)'"cnt)) ; 
end 

pbm= (pb''m)*pbm; 

%(*  calculate  the  value  for  P^  for  d  even,  equqation  (36)  *) 
pd=(2*sqrt(pbm*(l-pbm))),"'d; 

%(*  calculate  P^,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,pe=l;,end 
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i 


%(*  calculate  the  integration  for  q(ri),  equation  (38)  *) 
inc=(l“pe)'^L*fr(l,o)*fr(2,a)*:fr(2,b)*fr(2,c)*fr(2,f); 
inc=  mc*fr(2,g)*fr(2,h)*fr(2,k);. 
sum=sum+inc; 
end 
end 
end 
end 
end 
end 
end 
end 

sura=sum*delta(2)^(nuser-l)*delta(l); 

end 


The  source  code  for  the  numerical  results  of  Table  IV  is  listed  below.  For  ring 
two,  Rj  varies  from  25  km  to  38.5  km,  and  R2  varies  from  0  km  to  50  km.  This 
program  is  used  to  calculate  q(8).  To  calculate  q(2),  q(3),...,  q(n),  the  number  of 
users  (nuser)  is  equal  to  2,  3,...,  n,  respectively. 
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%(*  this  program  will  calculate  q(8)  for  ring  two  *)  , 

%(*  the  . coefficient  values  for  d  and  *) 
d  =[10  12  14  16  18  10]; 
ad=[ll  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  ’') 
interv=5; 

%(*  the  values  for  R,  and  R,  *) 
rl=[25  0]; 
r2=[38.5  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 
alpha =4; 

L=400; 
nuser=8; 
ml=[l  2  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2) = (r2(2)-rl(2))/interv; 

r(l,l)=rl(l)+delta(lV2; 

r(2,l)=rl(2)+delta(2)/2; 

for  i=2:interv 

r(l,i) = r(l,i-l) + delta(l); 
end 

for  i=2:interv 

r(2,i) =r(2,i-l) + delta(2); 
end 

%(*  calculate  the  probability  density  function  f^  equation  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(l)"2-rl(l)"2); 

end 

for  i=l:interv 

fr(2.i)=2*r(2.i)/(r2(l)^-rl(2)"2); 

end 
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%(*  this  loop  calculates  q(8)  for  all 
%value  of  M=  1,  2,  and  4,  equation  (38)  *) 
for  i2=l;3 
in=ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:interv 
for  h=l:interv 
for  k=l;interv 
for  o  =  l:interv 

index=[a  b  c  f  g  h  k  o]; 
muy=0; 

%{*  calculate  value  for  m,  equation  (27)  *) 
for  il=l:nuser-l 

muy=muy+ (r(l,index(l)/r(2,index(il  +  l)))'^alpha; 
end 

muy = 1/sqr  t(l + muy); 

%{*  calculate  value  for  K  equation  (25)  *) 
pb=(l-muy)/2 

%(*  calculate  the  value  for  equation  (26)  *) 

pbm=0; 

for  cnt=0:ni-l 

pbm=pbm+  (comb((m-l +cnt),cnt)*(((l+niuy)/2)'^cnt)); 
end 

pbm = (pb^m)  *pbm; 

%(*  calculate  the  value  for  for  d  even,  equqation  (36)  *) 
pd = (2*sqrt(pbm*(l-pbm)))."'d; 

%(*  calculate  P,.,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,pe=l;,end 
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%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc=(l-pe)^L*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f); 
inc = inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum=  sura + inc; 
end 
end 
end 
end 
end 
end 
end 
end 

sura=sura*delta(2)^(nuser-l)*delta(l); 

end 


The  source  code  for  the  numerical  results  of  Table  V  is  listed  below.  For  ring 
three,  varies  from  38.5  km  to  50  km,  and  R2  varies  from  0  km  to  50  km.  Tliis 
program  is  used  to  calculate  q(8).  To  calculate  q(2),  q(3),...,  q(n),  the  number  of 
users  (nuser)  is  equal  to  2,  3,...,  n,  respectively. 


%(*  this  program  will  calculate  q(8)  for  ring  three  *) 

%(*  the  coefficient  values  for  d  and  a^  *) 
d  =[10  12  14  16  18  10]; 
ad=[ll  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  and  R2  *) 
rl= [38.5  0]; 
r2=[50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 
alpha =4; 

L=400; 
nuser=8; 
ml=[12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l) =(r2(l)-rl(l))/interv; 

delta(2) = (r2(2)-rl(2))/interv; 

r(l,l)=rl(l)+delta(l)/2; 

r(2,l)=rl(2)+delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)+delta(l); 

end 

for  i=2:interv 

r(2,i) =r(2,i-l) + delta(2); 
end 

%(*  calculate  the  probability  density  function  fj^  eq  '  ion  (7)  *) 
for  i=l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(ir2); 

end 

for  i=l:interv 

fr(2,i)=2*r(2,i)/(r2(l)"2-rl(2)"2); 

end 


%(*  this  loop  calculates  q(8)  for  all 
%value  of  M=l,  2,  and  4,  equation  (38)  *) 
for  i2=l:3 
ni=ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l;interv 
for  f=l:mterv 
for  g=l:interv 
for  h=l:interv 
for  k=l;interv 
for  o=l:interv 

index=[a  b  c  f  g  h  k  o]; 
muy=0; 

%(*  calculate  value  for  equation  (27)  *) 
for  il  =  l:nuser-l 

muy =muy+  (r(l,index(l)/r(2,index(il  +  l)))'^alpha; 
end 

muy = l/sqrt(l  +muy); 

%(*  calculate  value  for  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  equation  (26)  *) 

pbm=0; 

for  cnt=0:m-l 

pbm=pbm+  (comb((m-l +cnt),cnt)*(((l+muy)/2)'^cnt)); 
end 

pbm=(pb''m)*pbm; 

%(*  calculate  the  value  for  P^  for  d  even,  equqation  (36)  *) 
pd = (2*sqrt(pbm*(l-pbm)))^d; 

%(*  calculate  P^,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,pe=l;,end 
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%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc=(l-pe)^L*fr(l,o)*fr(2,a)*£r(2,b)*fr(2,c)*fr(2,f); 
inc = mc*fr(2,g)*fr(2,h)*fr(2,k); 
sum=sum+inc; 
end 
end 
end 
end 
end 
end 
end 
end 

sum = sum*  delta(2)^(nuser-l)  *delta(l); 
end 
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APPENDIX  B 


Equation  (36)  is  used  to  calculate  q(n)  since  computing  q(n)  for  n  >  6  from 
equation  (30)  is  prohibitively  time-consuming.  The  following  programs  use  (30)  for 
Pj  to  calculate  q(n).  The  following  programs  are  listed  for  reference  purposes.  The 
source  codes  and  the  results  are  listed  below.  Since  q(n)  for  n  >  6  takes  a  long  time 
to  calculate,  the  program  is  written  to  calculate  only  up  to  q(6). 

The  number  of  users  (nuser)  is  initialized  to  2,  3 .  6  to  calculate  q(2),  q(3), 

...,  q(6),  respectively.  In  the  listed  program  the  number  of  users  (nuser)  is  equal  to 
six.  For  a  one-ring  network,  Rj  and  R2  vary  from  0km  to  50km.  The  source  code  and 
its  simulation  results  are  listed  below  and  in  Table  VI,  respectively. 
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%(*  this  program  calculates  q(6)  for  one  ring  netv/ork  *) 

%(•  the  coefficient  values  for  d  and  a<j  *) 
d  =[10  12  14  16  18  10]; 
ad=[ll  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Ri  and  R2  *) 
rl=[0  0]; 
i2=[50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *). 
alpha =4; 

L=400; 
nuser=8; 
ml=[l  2  4]; 

%(*  calculate  the  grid  points  of  integration  •) 

delta(l) = (r2(l)-rl(l))/interv; 

delta(2) = (r2(2)-rl(2))/interv; 

r(l,l)=rl(l)+delta(l)/2; 

r(2,l)=rl(2)+delta(2)/2; 

for  i=2:interv 

r(l,i) =r(l,i-l) + delta(l); 
end 

for  i=2:interv 

r(2,i)=r(2,i-l)+delta(2); 

end 

%(*  calculate  the  probability  density  ftinction  ^  equation  (7)  *) 
for  i=l:interv 

fr(l,i)=2*r(l,i)/(r2(l)''2.rl(ir2); 

end 

for  i=l;interv 

fr(2,i)=2*r(2,i)/(r2(l)"2-rl(2)"2); 

end 
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%(*  tliis  loop  calculate  q(6)  for  all 
%value  of  M=l,  2,  and  4,  equation  (38)  *) 
for  i2=l;3 
m=inl(i2); 
sum=0; 
for  a=l:interv 
for  b=l;interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:interv 
for  o=l:interv 

index = [a  b  c  f  g  o]; 
muy=0; 

%(*  calculate  value  for  n,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy = muy  +  (r(l,index(l)/r(2,index(i  1 + l)))^alpha; 
end 

inuy=l/sqrt(l+muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb=(l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm=0; 

for  cnt=0:ni-l 

pbm  =pbm+ (comb((m-l +cnt),cnt)*(((l +muy)/2)''cnt)); 

end 

pbni=(pb^m)*pbm; 

%{*  calculate  the  value  for  Pj  for  d  even,  equqation  (30)  *) 
for  q=l:6 

pd(q) = 0.5*comb(d(q),d(q)/2)*pbm^(d(q)/2); 
pd(q) =pd(q)*(l-pbm)^(d(q)/2); 

for  e = 1 + d(q)/2:d(q) 

pd(q) = pd(q) + (comb(d(q),e)  *  (pbrn^'e)  *  (1-pbm)  (d(q)-e)) ; 
end 
end 


%(*  calculate  P^,  equation  (29)  *). 

pe=ad*pd’; 

if  pe>li  pe=l;,  end 


%(*  calculate  the  integration  for  q(ri)j  equation  (38)  *) 
inc=(l-pe)"L*fr(l,o)*ff(2,a)*fr(2,b)*fr(2;c)*fr(2,f); 
inc=inc*fr(2,g); 
sum=sum+inc; 
end 
end 
end 
end 
end 
end 

sum=sum*delta(2)'^(nuser-l)*delta(l); 

end 


Table  VI.  SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  ONE-RING  NETWORK 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

M=1 

1 

.228 

.138 

.096 

.075 

.059 

M=2 

1 

.412 

.238 

.159 

.123 

.096 

M=4 

1 

.636 

.442 

.314 

.241 

.197 

The  source  code  and  the  simulation  results  are  listed  below  and  in  Table  VII, 
respectively.  For  ring  one  Rj  is  varied  from  0  km  to  25  km,  and  R2  is  varied  from  0 
km  to  50  km.  This  source  code  is  used  to  calculate  q(6).  To  compute  q(2), 
q(3),...,q(n),  the  number  of  users  (nuser)  is  2,3,...,n  respectively. 
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%{*  this  program  calculates  q(6)  for  ring  one  *) 

%(*  the  coefficient  values  for  d  and  a^  *) 
d  =[10  12  14  16  18  10]; 
ad =[11  38  193  1331  7275  40406]; 

%(*  the  number  of  inteival  for  integration  *) 
interv=5; 

%(*  the  values  for  and  R2  *) 
rl  =  [0  0]; 
r2=[25  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  .ml  for  M 
alpha =4; 

L=400; 
nuser=8; 
ml=[12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l) = (r2(l)-r  1  ( 1)) /interv; 

delta(2) = (r2(2)-rl(2))/interv; 

r(l,l)=rl(l)+delta(l)/2; 

r(2,l)=rl(2)+delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)+delta(l); 

end 

for  i=2:interv 

r(2,i)  =r(2,i-l) + delta(2); 
end 

%(*  calculate  the  probability  density  function  f^  equation  (7)  •) 
for  i=l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(ir2); 

end 

for  i=l;interv 

fr(2,i)=2*r(2,i)/(r2(ir2-rl(2r2); 

end 


%(*  this  loop  calculates  q(6)  for  all 
%value  of  M=l,  2,  and  4,  equation  (38)  *) 
for  12=1:3 
in=nil(i2); 
suin=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:iriterv 
for  o=l:interv 

index = [a  b  c  f  g  o]; 
muy=0; 

%(*  calculate  value  for  /i,  equation  (27)  *) 
for  il=l:nuser-l 

inuy =niuy + (r(l,index(l)/r(2,index(il + l)))'^alpha; 
end 

muy = l/sqrt(l+ muy); 

%(*  calculate  value  for  P^,  equation  (25)  *) 
pb=(l-muy)/2 

%(*  calculate  the  value  for  P,,jvi>  equation  (26)  *) 

pbm=0; 

for  cnt=0:ni-l 

pbm =pbm + (conib((m-l + cnt),cnt)*(((l +muy)/2)^cnt)); 
end 

pbm = (pb^m)*pbm; 

%(*  calculate  the  value  for  for  d  even,  equqation  (30)  *) 
for  q=l:6 

pd(q)=0.5*comb(d(q),d(q)/2)*pbm'"(d(q)/2); 
pd(q)  =pd(q)*(l-pbm)^(d(q)/2); 
for  e=H-d(q)/2:d(q) 

pd(q) =pd(q) + (comb(d(q),e)*(pbm^e)*(l-pbm)^(d(q)-e)); 
end 
end 
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%{*  calculate  P^.,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,  pe=l;,  end 

%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 

inc=(l-pe)X*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f); 

inc=inc*fr(2,g); 

.  jm=sum+inc; 
end 
end 
end 
end 
end 
end 

sum=sum*delta(2)^(nuser-l)*delta(l); 

end 


Table  VII.  SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  RING  ONE 


q(i) 

q(2) 

q(3) 

q(4) 

mm 

q(6) 

M  =  1 

1 

.751 

.505 

.355 

.2838 

.238 

M=2 

1 

.881 

.752 

.624 

.5169 

.437 

M=4 

1 

.916 

.836 

.761 

.688 

.617 

The  source  code  and  the  simulation  results  are  listed  below  and  in  Table  VIII, 
respectively.  For  ring  two,  Rj  is  varied  from  25  km  to  38.5  km,  and  R2  is  varied  from 
0  km  to  50  km.  This  program  is  used  to  calculate  q(6).  To  calculate  q(2),  q(3),..., 
q(n),  the  number  of  users  (nuser)  is  2,  3,...,  n  respectively. 
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%{*  thii,  program  will  calculate  q(6)  for  ring  two  *) 
the  coefficient  values  for  d  and  a^  *) 
d  =[10  12  14  16  18  10]; 
ad = [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *); 

interv=5;  '  ^ 

%{*  the  values  for  R,  and  R,  *) 
rl=[25  0]; 
r2=[38.5  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 
alpha =4; 

L=400; 
nuser=8; 
ml=[12  4]; 

%(*  calculate  the  grid  points  of  integration  •) 

delta(l) = (r2(l)-rl(l))/interv; 

delta(2) = (r2(2)-rl(2))/interv; 

r(l,l)=rl(l)+delta(l)/2; 

r(2,l)=rl(2)+delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)+delta(l); 

end 

for  i=2:interv 

r(2,i)=r(2,i-l)+delta(2); 

end 

%(*  calculate  the  probability  density  function  f^  equation  (7)  *) 
for  i=l;interv 

fr(l,i)=2’r(l.i)/(r2(ir2-rl(ir2); 

end 

for  i=l:interv 

fr(2.i)=2*r(2.i)/(r2(l)"2-rl(2)"2); 

end 


%(*  this  loop  calculates  q(6)  for  all 
%value  of  M=l,  2,  and  4,  equation  (38)  *) 
for  12=1:3 
m=ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l;interv 
for  g=l:interv 
for  o=l:interv 

index = [a  b  c  f  g  o]; 
muy=0; 

%(*  calculate  value  for  n,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy = muy + (r(l,index(l)/r(2,index(i  1 + l)))'^alpha; 
end 

muy = l/sqrt(l + muy); 

%(*  calculate  value  for  P^,  equation  (25)  *) 
pb=(l-muy)/2 

%(*  calculate  the  value  for  equation  (26)  *) 

pbm=0; 

for  cnt=0:m-l 

pbm =pbm + (comb((m-l + cnt),cnt)*(((l +muy)/2)^cnt)); 
end 

pbm = (pb''m)’pbm; 

%(*  calculate  the  value  for  Pj  for  d  even,  equqation  (30)  *) 
for  q=l:6 

pd(q) = 0.5*comb(d(q),d(q)/2)*pbm'"(d(q)/2); 
pd(q) =pd(q)*(l-pbm)^(d(q)/2); 
for  e=  1+ d(q)/2:d(q) 

pd(q) =pd(q) + (comb(d(q),e)*(pbm^e)*(l-pbm)'"(d(q)-e)); 
end 
end 
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%(*  calculate  Pg,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,  pe=l;,  end 


%(*  calculate  the  integration  for  q(n),  equation.  (38)  *) 
inc=(l-pe)"L’fr(l,o)*£r(2,a)*fr(2,b)*fr(2,c)*fr(2,f); 
inc=inc*fr(2,g); 
sum=sun+inc; 
end 
end 
end 
end 
end 
end 

sum=sum*delta(2)^(nuser-l)*delta(l); 

end 


Table  VIII.  SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  RING  TWO 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

II 

1 

.055 

.0018 

0 

0 

0 

M=2 

1 

.33 

.1021 

.0227 

.0052 

.0011 

M=4 

1 

.487 

.2971 

.1534 

.065 

.0293 

The  source  code  and  the  simulation  results  are  listed  below  and  in  Table  IX, 
respectively.  For  ring  three,  Rj  is  varied  from  38.5  km  to  50  km,  and  R2  is  varied 
from  0  km  to  50  km.  This  program  is  used  to  calculate  q(6).  To  calculate  q(2), 
q(3),...,  q(n),  the  number  of  users  (nuser)  is  2,  3,...,  n  respectively. 


%(*  this  program  will  calculate  q(6)  for  ring  three  *) 

%(*  the  coefficient  values  for  d  and  a<j  *) 
d  =[10  12  14  16  18  10]; 
ad =[11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  and  Rj  *) 
rl  =  [38.5  0'j; 
r2=[50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 
alpha =4; 

L=400; 
nuser=8; 
ml=[12  4]; 

%(*  calculate  the  grid:  points  of  integration  *) 
delta(l) = (r2(l)-rl(l))/interv; 
deita(2) =(r2{2)-rl(2))/interv; 
r(l,l)=rl(l)+delta(l)/2;  ' 
r(24)=rl(2)+delta(2)/2; 
for  i=2:interv 

r(l,i) =r(l,i-l) + delta(l); 
end 

for  i=2:interv 

r(2,i) = r(2,i-l) + delta(2); 
end 

%(*  calculate  the  probability  density  function  f^  equation  (7)  *) 
for  i=l:interv 

fr(l,i)=2*r(l,i)/(r2(l)"2-rl(l)"2); 

end 

for  i=l:interv 

fr(23i)=2*r(2.i)/(r2(l)''2-rl(2)'^); 

end 


%(*  tliis  loop  calculates  q(6)  for  all 
%value  of  M=l,  2,  arid  4,  equation  (38)  *) 
for  i2=l:3 
m=ml{i2); 
sum=0; 
for  a=l:interv 
for  b-l:interv 
for  c=l;mterv 
for  f=l:interv 
for  g=r:mterv 
for  o=l:interv 

index=[a  b  c  f  g  o]; 
muy=0; 

%(*  calculate  value  for  n,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy = muy + (r(l,index(l)/  r(2,index(i  1  +  l)))'^alpha; 
end 

muy = l/sqrt(l + muy); 

%(*  calculate  value  for  Pj,,  equation  (25)  *) 
pb=(l-muy)/2 

%(*  calculate  the  value  for  equation  (26)  *) 

pbm=0; 

for  cnt=0:m-l 

pbm =pbm+  (comb((m-l + cnt),cnt)*(((l +muy)/2)'"cnt)); 
end 

pbm  =  (pb''m)*pbm; 

%{*  calculate  the  value  for  P^  for  d  even,  eqiiqation  (30)  *) 
for  q= 1:6 

pd(q) = 0.5*comb(d(q),d(q)/2)’pbm'"(d(q)/2); 
pd(q)  =pd(q)*(l-pbm)^(d(q)/2); 
for  e = 1 + d(q)/2:d(q) 

pd(q) =pd(q) + (comb(d(q),e)*(pbm^e)*(l-pbm)^(d(q)-e)); 
end 
end 
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%(*  calculate  P^,  equation  (29)  *) 

pe=ad*pd’; 

if  pe>l,  pe=l;,  end 

%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc = (l-pe)X*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,Q; 
incline’*  fr(2,g); 
sum = sum + inc; 
end 
end 
end 
end 
end 
end 

sum=sum*delta(2)^(nuser-l)*delta(l); 

end 


Table  IX.  SIMULATION  RESULTS  FOR  THE  PROBABILITY  OF 
CORRECT  PACKET  CAPTURE  FOR  RING  THREE 


q(i) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

M  =  1 

1 

0 

0 

0 

0 

0 

M=2 

1 

.0474 

0 

0 

0 

0 

li 

1 

.2578 

.0447 

.003 

.002 

0 

Plots  of  q(n)  versus  n  for  one-ring  and  three-ring  systems  are  shown  in  Figure 
11,  12,  13,  and  14. 
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Figure  11.  Probability  Of  Correct  Packet  Capture  For  One-Ring  Network 
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Total  Number  of  P 


Figure  12.  Probability  Of  Correct  Packet  Capture  For  Ring  One 
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Figure  13.  Probability  Of  Correct  Packet  Capture  For  Ring  Two 
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Figure  14.  Probability  Of  Correct  Packet  Capture  For  Ring  Three 
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APPENDIX  C 


%(*  initialization  *) 
nuser=15; 

p0=0.1; 

pr= 0.002; 

q=[0  1  0.9042  0.8099  0.7197  0.6350  0.5580  0.4901  0.4615  0.4285]; 

%(*  steady  state  value  of  q(n)  *) 
for  i=ll;nuser+4 
q(i)= 0.3852; 
end 

%(*  set  N=nuser,  calculate  it  is  called  pi(n,m)  *) 

N  =nusef; 
for  n=0:N 
nl=n+l; 
for  m=0:N 
ml=m+l; 

if  m<n-l,pi(nl,ml)=0; 

elseif  m  =  =  n-l 
suml=0; 
for  i=0;n 
il=i+l; 

sunil=suml+comb(n,i)*pr''i*(l-pr)^(n-i)*q(il); 

end 

pi(nl,ml) = (l-p0)'^(N-n)*suml; 

elseif  m>  =n 
suml=0; 

templ=comb(N-n,N-m)*(l-p0)''(N-m-l)*p0"'(m-n); 
for  i=0:n 
il=i+l 

term = comb(n,i)  *  (l-pr)""(n-i)  *pr; 

term=term*((l-p0)*(l-q(il + m-n)) + (N-m)*p0*q(il +m-n+  l)/(m-n+ 1)); 
suml=suml+term; 
end 

pi(nl,ml)=templ*suml; 


%(*  calculate  7r„*,  it  is  called  pis(n)  *) 
for  n=l:N 
nl=n+l; 
suml=0; 
for  i=0:n-l 
il=i+l; 

suml =suml  +  (pis(il)*pi(il,nl-l)); 
end 

temp=pis(nl-l)-suml; 
pis(nl) = (l/pi(nl,nl-l))  *  temp; 
if  pis(nl)<0,pis(nl)=0;,end 
end 

%(*  calculate  7r„,  it  is  called  pin  *) 

den=sum(pis); 

pin=pis/den; 

%(*  calculate  fn  *) 
for  n=0:N 
nl=n+l; 
temp=0; 
total =0; 
for  k=0:N-n 
kl=k+l; 
suml=0; 
for  i=0:n 

suml = suml  +  (comb(n,i)*pr‘^i*(l-pr)'^(n-i)*q(i + kl); 
end 

temp = comb  (N-n,N-n-k)  *p0""k(  1  -pO)''(N-n-k) ; 
term = suml*  temp; 
total = total + term; 
end 

f(nl)= total; 
end 

%(*  calculate  s  *) 
suml=0; 
fori=0:N 
il=i+l; 

suml=suml+(pin(il)*f(il)); 

s(il)=suml; 

end _ 


%{*  calculate  the  Comb(n,m)  *) 
function  y=corab(n,m) 
y =xfact(n+  l)/(xfact(m+  l)*xfact(n-m+ !)); 

%(*  calculate  the  factorial,  it  is  call  fact(n)  *) 
fiinction  y=fact(x) 
ifx==0 
y=l; 
else 

y=i; 

for  i=l:x 
y=y*i; 
end 
end 

%(’  the  variable  xfact  *) 

x=0:52 

for  i  =  l:53 

xfact(i)  =fact(x(i)); 
end 

global  xfact 
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